Optimal working points for continuous-variable quantum channels 
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The most important ability of a quantum channel is to preserve the quantum properties of trans- 
mitted quantum states. We experimentally demonstrate a continuous- variable system for efficient 
benchmarking of quantum channels. We probe the tested quantum channels for a wide range of 
experimental parameters such as amplitude, phase noise and channel lengths up to 40 km. The data 
is analyzed using the framework of effective entanglement. We subsequently are able to deduce an 
optimal point of operation for each quantum channel with respect to the rate of distributed entan- 
glement. This procedure is a promising candidate for benchmarking quantum nodes and individual 
links in large quantum networks of different physical implementations. 
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It is envisaged that a quantum network would consist 
of many nodes and channels for processing and distribu- 
tion of information [1]. Such a network requires the abil- 
ity to perform well-controlled operations on single quan- 
tum systems and to generate entanglement between sev- 
eral quantum systems. In the past, considerable efforts 
were made to implement, characterize and optimize such 
local operations [2]. As the complexity of single nodes 
increases, another crucial aspect to consider is the trans- 
fer of quantum information within and between them. In 
quantum networks, as well as quantum computers, char- 
acterization of the individual nodes and channels via pro- 
cess tomography [3] or a direct benchmarking procedure 
is essential to get a performance estimate of the whole 
system [4, 5]. Moreover, if we can analyze nodes and 
channels quantitatively, then, by tuning parameters, we 
can optimize their operational working points. 

A common way to evaluate the performance of quan- 
tum devices is to compare them against analogous clas- 
sical devices using fidelity-based benchmarks [6], even 
though this approach can be experimentally costly. 
While fidelity may allow us to distinguish devices from 
their classical counterparts, it is not clear how to use 
it to quantitatively characterize devices. Results from 
the field of quantum key distribution (QKD) can pro- 
vide another possible way to quantitatively character- 
ize quantum networks. In QKD, one aims to generate 
a secret key from quantum-correlated data by measur- 
ing quantum states shared between distant parties. The 
quantum states are exchanged over a quantum channel. 
This channel is an important component and requires 
great care when being established and characterized [7]. 
However, there are other less demanding applications for 
quantum channels, and the specific requirements for the 
channel will depend on our needs. This allows us to iden- 



tify some minimal requirements [8] , based on the preser- 
vation of entanglement, for a channel to be useful for 
quantum communication. 

From a more general perspective, the distribution of 
quantum correlations in quantum networks can be effi- 
ciently investigated using the framework of effective en- 
tanglement [4, 5, 9]. Although our approach applies gen- 
erally, we focus here on continuous- variable systems, such 
as systems based on quantum optics, collective spin ex- 
citations [10], surface plasmons [11], optomechanics [12] 
and cavity quantum electrodynamics [13]. In contrast to 
fidelity-based benchmarks, the framework of effective en- 
tanglement aims to use minimal experimental resources 
while providing a quantitative characterization of quan- 
tum devices. So far, experiments have been performed to 
witness effective entanglement in a fiber [14] and a free- 
space [15] channel, as well as a benchmark performed on 
a quantum memory [16]. It is now desireable to have 
a benchmark system that is capable of covering a wide 
range of parameters in order to find the optimal point of 
operation for each quantum device. 

In this Letter, we demonstrate an experimental imple- 
mentation of a continuous-variable system for efficient 
benchmarking of quantum channels and devices. This 
benchmark system allows us to identify the optimal point 
of operation of a given quantum channel with respect to 
the rate of transmitted entanglement. We find the opti- 
mum by varying the quantum state amplitude and mea- 
suring the excess noise for fiber channels with a length 
of up to 40 km. The entanglement rates are calculated 
using the framework of effective entanglement. Further- 
more, we demonstrate the impact of phase noise on the 
ability to distribute entanglement, showing that after a 
certain amount of dephasing the detectable quantum cor- 
relations drop to zero. 
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The primary goal of quantum benchmarking is to op- 
erationally differentiate a given channel/device from the 
class of measure and prepare (MP) channels. MP chan- 
nels involve first a measurement, followed by classical 
storage/transmission of the measurement result, and fi- 
nally the preparation of output states based on this clas- 
sical data. Due to the classical step, MP channels are not 
sufficient for quantum communication protocols. Accord- 
ingly, any channel which cannot be distinguished from 
MP operation will not be usable in a quantum commu- 
nication context. 

The set of MP channels is mathematically equivalent to 
the set of entanglement-breaking channels [17]. Thus, if 
we show that a channel can preserve entanglement, then 
that channel has passed our quantum benchmark. In 
fact, by using a source-replacement (or prepare and mea- 
sure) description, the entanglement-preserving capabili- 
ties of a channel can be efficiently tested without needing 
to use real entangled states. In such a scheme, we probe 
the channel with different test states (for the experiment, 
we use two coherent states |±a)). Theoretically, we imag- 
ine that these test states are part of a larger entangled 
state, given by [9, 18] 
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We represent the tested channel mathematically by a 
quantum map A, mapping from the input system A' 
to some output system B. To test the entanglement- 
preserving capabilities of the channel, we need to deter- 
mine whether or not the bipartite output state p jfe = 
(id ® A) [|^)(^| AA ,] is entangled. However, since the en- 
tangled state is a virtual construction, we do not have 
access to tomographically-complete information of this 
effectively entangled state. Our information is limited to 
only the reduced density matrix pA = Tr# (//£#), along 
with any measurements made on the physical output 
states Po/r ^he recm ced density matrix 
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is parameterized by the overlap s = (a\ —a) and remains 
constant during the protocol. Previously, we have used 
Stokes measurements [14] because they enable bench- 
marking under the adversarial scenario common in QKD. 
Here, we are interested in quantifying channels under less 
restrictive conditions. In this case, the conjugate quadra- 
ture operators x and p suffice (see [18] for discussion 
about Stokes vs quadrature measurements). Through 
measurement, we determine the mean value (z) and vari- 
ance Var(i) = (i 2 ) — (z) 2 of z G {x,p} for each output 
state. 

Using the overlap value and the obtained measure- 
ment results, the task is to determine whether p°£g is 
entangled. One method for accomplishing this is with 
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FIG. 1: (a) Quantum benchmarking system using a prepare- 
and-measure-type setup. Effective entanglement can be quan- 
tified using Bob's measurement data of the signal states, (b) 
Schematic drawing of the experimental setup. The abbrevi- 
ated components are beamsplitters (BS), the Mach-Zehnder 
modulator (mod.), an attenuator (att.), an optical isolator 
(iso.), a polarizing beamsplitter (PBS), quarter- wave plate 
(QWP) and half- wave plate (HWP). 



a witness-like procedure [9, 18]. Alternatively, we can 
minimize an appropriate entanglement measure over all 
bipartite states consistent with the available information 
[5, 19]. In addition to classifying entangled states [5], 
the entanglement measure approach gives us the abil- 
ity to also quantitatively describe and compare different 
channels. Although we use the concept of effective entan- 
glement to simplify our benchmarking scheme, we note 
that it is entirely representative of what one would find 
had the same data been generated from a real entangled 
state. 

Our entanglement measure cannot be chosen arbitrar- 
ily; the quantification procedure outlined in [5, 19] re- 
quires that the measure satisfies some mild monotonic- 
ity properties. In addition, the measure should ideally 
be efficiently computable, for instance via a semidefinite 
program. In this paper, we will use the negativity [20- 
22], which fulfills the required conditions. It will also be 
convenient to convert our entanglement content to the 
logarithmic negativity [23], another measure which is in 
one-to-one correspondence with the negativity. The log- 
negativity has the useful property of being additive, al- 
lowing us to quantify entanglement built up over repeated 
channel use. It also has an operational interpretation in 
terms of the entanglement cost under a certain class of 
operations [24]. 

We implemented the benchmark using an experimen- 
tal setup similar to that demonstrated in [14], consisting 
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of three parts: a sender (Alice), a receiver (Bob) and the 
channel to be tested (which is simply "plugged in" to 
the benchmark system) (Fig. 1). Alice and the quantum 
channel are fiber-based, whereas Bob uses free space op- 
tics. We use a 1550 nm laser as a light source, operated 
at a repetition rate of 1 MHz, producing pulses of 20ns. 
The pulses are then split up asymmetrically (BS1) into 
a local oscillator (LO) and a signal line. The larger por- 
tion of the beam is directed to the LO line to provide 
a strong phase reference for the weak signal beam. The 
signal preparation line is composed of a delay fiber, a 
Mach-Zehnder modulator (MZM), an attenuator and a 
polarization controller. The delay fiber causes a 0.5 (is 
delay between signal and LO pulses to prevent scatter- 
ing of photons from the LO into the signal mode during 
their propagation through the channel. The MZM gener- 
ates coherent states with opposite phases and ir while 
the attenuator reduces the state amplitude to a quan- 
tum level. A polarization controller (PCI) prepares the 
signal states' polarization orthogonal to the LO pulses. 
After recombining the two modes with a beamsplitter 
(BS2) another polarization controller (PC2) is adjusted 
such that LO and signal can be split up again on Bob's 
side. The two modes are then sent through an isolator 
to the quantum channel, after which they enter Bob's 
setup. There, the time multiplexing delay is reversed to 
allow for temporal mode matching of signal and LO. The 
signal states are then characterized using double homo- 
dyne detection. Additionally, the signal pulses in Alice 
and the LO in Bob are monitored. 

We employ a binary alphabet to probe the quantum 
channel. The alphabet consists of two weak coherent 
states: \a) and \—ct) ~ 0.5). This alphabet is cho- 
sen, because it is easy to implement (single MZM) and to 
post-process. After the states have passed through the 
quantum channel they enter Bob's receiver. We mea- 
sured a receiver efficiency J] Bob of ~ 73%, which consists 
of the optical losses in Bob's setup, the homodyne de- 
tector efficiencies and the interferometric visibility be- 
tween the signal and LO mode. Note that the total 
transmission T used later on includes the receiver effi- 
ciency TjBob and the losses in the quantum channel. Bob's 
double homodyne detection measures the joint probabil- 
ity distribution of the conjugate quadrature operators 
x and p, also known as the Q-function [25]. We ob- 
tain the marginal distribution of the Q-function for each 
quadrature to calculate its first and second moments. 
From the first moments we compute the complex num- 
ber j3 = (a) 1 for which we use the symmetric definition 
(a) = + i(p)). Knowing the total transmission T 

we are able to calculate the overlap (a\—a) of the co- 




1 a denotes the quantum mechanical annihilation operator for the 
light mode under investigation. 
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FIG. 2: Quantification of effective entanglement for a 20 km 
(top) and a 40 km (bottom) fiber channel. The negativity 
values in the shaded region were calculated by assuming all 
variances equal [5]. Since the experimental variances have 
small differences in x and p of the alphabet, the theoretical 
negativity values here serve only as a reference. We observe 
that the 20 km channel can tolerate more excess noise whereas 
the 40 km channel shows to be at its limit. Note that the dis- 
played negativity is a lower bound to the actual transmitted 
entanglement and that the statistical error is contained within 
the data points. 



herent states on Alice's side. The second moments of 
the Q-function marginals, minus a constant, correspond 
to the quadrature variances Var(x) and Var(p). We nor- 
malize the measured variances of our signal states to that 
of the coherent vacuum state, where we chose the shot 
noise level to be equal to 1. Any additional variance is 
considered to be excess noise added by the channel. 

The parameters we varied are the transmission of our 
quantum channel, the sending amplitude \a\ of the co- 
herent test states and the magnitude of the artificially 
induced phase noise. We performed the benchmark for 
two standard fibers with lengths of 20 km (T « 24%) 
and 40 km (T « 9%) serving as the quantum channel. 
We probed the channels for different sending amplitudes 
ranging approximately from < \a\ < 1 and found a 
nearly constant quadrature variance of « 1.02 (20 km 
fiber) and ~ 1.04 (40 km fiber). We were able to wit- 
ness and quantify the preservation of entanglement for 
both quantum channels for certain amplitudes (Fig. 2). 
Using a sending rate of 0.875 MHz and the quantified 
log- negativity we calculated the rate of distributed entan- 
glement for both channels (Fig. 3). We find maximum 
rates of 166,000 log-neg units/s for a 20 km and 15,000 
log-neg units/s for a 40 km channel. For comparison, a 
maximally-entangled two-qubit state has 1 log-neg unit 
of entanglement. We were thus able to determine the 
optimal point of operation for the two quantum channels 
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FIG. 3: Minimal rate of distributed entanglement for a 20 
km (top) and a 40 km (bottom) channel, obtained using ex- 
perimentally measured variances. Temporal variability in the 
entanglement transmission rates are determined by including 
la, 2ct and 3cr error bars in the optimization. In the exper- 
iment the temporal variability corresponds to fluctuations in 
the channel. 



TABLE I: Deduced optimal point of operation for the tested 
fiber quantum channels. 





T 


\a\ 


~ Var(x,p) 


negativity log-neg rate 


20 km 


24% 


0.3 


1.02 


0.07 0.19 166,000 


40 km 


9% 


0.17 


1.04 


0.006 0.017 15,000 



(Table I). 

Furthermore, we investigated the impact of excess 
noise on the channel and the effective entanglement 
benchmarking protocol. Here, our phase calibration pro- 
cedure provides access to a tunable noise parameter. As 
no phase lock is present between signal and LO pulses, 
their relative phase drifts freely whenever they propagate 
through different spatial paths of the setup. Without any 
post-processing the signal states thus arrive effectively 
phase-randomized. However, we experimentally find that 
their relative phase remains fairly constant during a 250 
lis time frame. Alice therefore sends along bright cali- 
bration pulses to provide a classical phase reference for 
Bob. These calibration pulses are used to rotate the sig- 
nal pulses in the measured phase space to match Alice's 
preparation frame of reference. This procedure effectively 
reverts the phase-randomization [26]. 

The phase noise may now be tuned, by gradually weak- 
ening the calibration pulses and thus diminishing the 
phase information available for the phase space rotation. 
We therefore expect a gradual phase-diffusion of the sig- 
nal states with decreasing calibration pulse amplitude. 
This phase-diffusion is reflected by an increased quadra- 



FIG. 4: Examination of entanglement throughput versus 
phase noise (the phase noise tuning parameter is defined as 
6 — |a ca i|)- The transition between quantum and classical do- 
mains can be observed for increasing phase noise (diminishing 
calibration amplitude). This is reflected in increased quadra- 
ture variances (top) and subsequently decreasing transmission 
rates (bottom). The rates have been obtained by comput- 
ing the log- negativity for the measured asymmetric variances. 
The channel was probed for \a\ « 0.5. (Inset figure: Illustra- 
tion of a phase- diffused coherent state.) 

ture variance of x and p. Below a certain calibration 
pulse amplitude, the signal states will be phase-diffused 
to such a degree, that the negativity drops to zero and 
entanglement can no longer be witnessed. Fig. 4 shows 
this effect for a 20 km fiber channel. 

The theoretical curves for Fig. 4 were obtained by 
numerically calculating the variance of the Q-function 
marginals of phase-diffused coherent states. The phase- 
diffused Q-function Q p d was obtained by calculating the 
weighted integral of Q-functions of coherent states with 
the same amplitude r but different phase: 

Q P d(/?) = - /(</>) | (re^/?) | 2 (3) 

n J(f>=0 

The amount of phase-diffusion of Q p d is determined 
by the phase probability distribution /(</>) of a coherent 
state, given by /(</>) = J Q C oh(^, 0) r&r. The amplitude 
of Qcoh is equal to the amplitude of the calibration states 
Bob receives, thus corresponding to the amount of phase 
information available to him. The asymmetry between 
Var(x) and Var(p) is merely a matter of choosing a fixed 
frame of reference in phase space for the simulated phase- 
diffused state (modulation along x). 

In conclusion, we demonstrated a versatile continuous- 
variable system for efficient benchmarking of quantum 
communication channels and devices. This procedure en- 
ables us to investigate the quantum throughput of vari- 
ous quantum systems independent of their physical real- 
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ization. We can probe the channel over a range of tun- 
able parameters, such as input state amplitude and phase 
noise, allowing us to determine optimal working points. 
Using the framework of effective entanglement, we deter- 
mined these points with respect to the rate of distributed 
entanglement. 
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